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We investigated the peculiarities of non-equilibrium charge states and spin configurations in the 
system of two strongly coupled quantum dots (QDs) weakly connected to the electrodes in the 
presence of Coulomb correlations. We analyzed the modification of non-equilibrium charge states 
and different spin configurations of the system in a wide range of applied bias voltage and revealed 
well pronounced ranges of system parameters where negative tunneling conductivity appears due to 
the Coulomb correlations. 
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I. INTRODUCTION 

Electron tunneling through the system of coupled 
quantum dots in the presence of strong Coulomb corre- 
lations seems to be one of the most interesting problems 
in the solid state physics. The present day experimental 
technique gives possibility to produce single QDs with a 
given set of parameters and to create coupled QDs with 
different spatial geometry fij, 0, @ , [H • Well known ver- 
tically aligned geometry [H , : 0] gives an opportunity to 
analyze non-stationary effects in various charge and spin 
configurations formation in the small size structures both 
theoretically [2l[ and experimentally Q. Lateral QDs 
are also actively studied experimentally during the last 
several years [J], but due to the technological problems 
they are mostly analyzed theoretically [5;] , [6] . Thereby 
the main effort in the physics of QDs is devoted to the 
investigation of non-equilibrium charge states and dif- 
ferent spin configurations due to the electrons tunneling 
; [III ; [HI , [H[ through the system of coupled QDs in 
the presence of strong Coulomb interaction. One of the 
most intensively studied problems in this field is tunnel- 
ing through the single QD , [TtJ , [3 and interacting 
QDs PJ], 0, [i^, [2(J, HU in the Kondo regime, which 
reveals rich physics for small bias voltage compared to 
the tunneling rates. Theoretical analysis of this problem 
usually deals with the Keldysh non-equilibrium Green- 
function formalism [22j.(23j. renormalization-group the- 
ory 12011 or specific approach suggested by Coleman 

m HI. HL m 

One of the most interesting results was obtained for the 
system of double QDs with on site Coulomb repulsion 
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in both of them [15|. The authors demonstrated that 
due to the presence of Coulomb correlations QDs can 
have a bistable behavior in the Kondo regime at zero 
bias voltage. 

Charge redistribution between different spin configura- 
tions in the system of two interacting QDs in the Kondo 
regime was regarded in |2l| . The authors considered 
the situation when the detuning between the energy lev- 
els in the QDs exceeds the dots coupling and on-site 
Coulomb repulsion is present only in a single dot. A 
new mechanism which leads to the transition from the 
singlet state in a weak coupling regime to a triplet state 
in a strong coupling regime was proposed. The authors 
demonstrated that interaction with continuous spectrum 
at zero bias modifies the energy of the singlet and triplet 
states and the situation when triplet state energy is lower 
than singlet one becomes possible. So careful analysis 
of tunneling processes through the system of interacting 
QDs in the Kondo regime reveals exciting physical phe- 
nomena. 

In the present paper we consider electron tunneling 
through the coupled QDs in the regime when applied bias 
can be tuned in a wide range and the on-site Coulomb 
repulsion can be comparable to the other system pa- 
rameters. We analyze all charge and spin configura- 
tions in the system of two strongly coupled quantum 
dots (QDs) weakly connected to the electrodes in the 
presence of Coulomb correlations in a wide range of ap- 
plied bias in terms of pseudo operators with constraint 

[I3,[H,d,0,[3,Eii For lar g e values of applied bias 
Kondo effect is not essential so we neglect any correla- 
tions between electron states in the QDs and in the leads. 
This approximation allows to describe correctly non equi- 
librium occupation of any single- and multi-electron state 
due to the tunneling processes. 

We revealed the presence of negative tunneling con- 
ductivity in certain ranges of the applied bias voltage 
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and analyzed the multiple charge redistribution between 
the two electron states with different spin configurations 
(singlet state and triplet state) as a function of the ap- 
plied bias voltage. 



II. THEORETICAL MODEL 

We consider a system of coupled QDs with the single 
particle levels £1 e 2 connected to the two leads. Such 
system can be described by means of two-impurities An- 
derson Hamiltonian where the impurities are the QD's 
24 1 , [2a| , [26| . The Hamiltonian can be written as: 



H = }^ af a a 1(7 ei + E ajg^^ + t/i"ig"i-<7 + 

a a 

+ U 2 n 2 an2-cr + ^T(a^ a a 2 a + a^ a a la ) (1) 



The operator creates an electron in the dot i with 
spin a, Si is the energy of the single electron level in 
the dot i and T is the inter-dot tunneling coupling, 
riic = OfoCLia and U\$ is the on-site Coulomb repul- 
sion of localized electrons. We'll consider for simplicity 
the situation when resonant tunneling between the QDs 
takes place and Coulomb repulsion is the same in the 
first and second QDs, consequently £ 1 = £2 = £0 and 
U\ = U 2 = U. Without the interaction with the leads 
all energies of single- and multi-electron states are well 
known: 

One electron in the system: two single electron states 
with energies e% = £0 ± T and wave function 



^ = — (|0t)|00)±|00>|0t» 



(2) 



Two electrons in the system: two states with the 
same spin a a and —a — a (triplet states has the spin 
projection Sz = ±1) with energies 2£ and four two- 
electron states with the opposite spins a — a and differ- 
ent configurations with energies Effj : 2e ; 2£ + U and 

Wave functions have the form: 



2£ ±^ 



^+4T 2 . 



rf° = a-|tt>|OO)-0-UO)|Ot> + 

+ 7-|ot>Uo} + Hoo>m) 



(3) 



Three electrons in the system: two three-electron 
states with energies E™fi = 3£o+£/±T and wave function 



Four electrons in the system: one four-electron state 
with energy Ejvi = 4eq + 2U and wave function 



i>i = \n)\n) 



(5) 



If coupled QDs are connected with the leads of the 
tunneling contact the number of electrons in the dots 
changes due to the tunneling processes. Transitions be- 
tween the states with different number of electrons in the 
two interacting QDs can be analyzed in terms of pseudo- 
particle operators with constraint on the physical states 
(the number of pseudo-particles) . Consequently, the elec- 
tron operator can be written in terms of pseudo- 
particle operators as: 



+ E V4,-^r CT + E ^™dj°~° + E tfVW 



(6) 



where f^ifa) and ip£(ip a )- are pseudo-fermion cre- 
ation(annihilation) operators for the electronic states 
with one and three electrons correspondingly. b + (b), 
d^(d a ) and tp + (tp)- are slave boson operators, which cor- 
respond to the states without any electrons, with two 
electrons or four electrons. Operators tj)^ n _ a - describe 
system configuration with two spin up electrons a and 
one spin down electron —a in the symmetric and asym- 
metric states. 

The constraint on the space of the possible system 
states have to be taken into account: 



n b 



■ E ™f ia + E "ST + E "V""o- + n v = 1 (7) 



Condition ([7]) means that the appearance of any two 
pseudo-particles in the system simultaneously is impos- 
sible. 

Electron filling numbers in the coupled QDs can be ex- 
pressed in the terms of the pseudo-particles filling num- 
bers: 



i 

+ E ^m-<7 + E 



E ; < 



13 y 



(8) 



Vw = • U o>u o> • (|o t>|oo> ± |oo>|o t» 

m = ±1 



Consequently, the Hamiltonian of the system can be 
(4) written in the terms of the pseudo-particle operators: 
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H = H +H tun (9) 

h = E e */-/- + E E™d+°°'df + 

+ E E Tn^ta^raa + E IVl ip+ip la + 
ma 

+ £ ka4 a c ka + J2(e pa - eV)c+ a c p(y 

ka per 
Htun = ^ T k( C ta C °- + C t c fccr) + ( fc ^ P) 



where e*, Effj , Effi and Eivi-are the energies of 
the single-, double-, triple- and quadri-electron states. 
£fc(p)o--is the energy of the conduction electrons in the 
states k and p correspondingly. c k ' p ^ )rT /ck(p)a are the cre- 
ation(annihilation) operators in the leads of the tunnel- 
ing contact. Tfc( p )-are the tunneling amplitudes, which 
we assume to be independent on momentum and spin. 
Indexes k(p) mean only that tunneling takes place from 
the system of coupled QDs to the conduction electrons 
in the states k and p correspondingly. 

Bilinear combinations of pseudo-particle operators are 



closely connected with the density matrix elements. So, 
similar expressions can be obtained from equations for 
the density matrix evolution but method based on the 
pseudo particle operators is more compact and conve- 
nient. The tunneling current through the proposed sys- 
tem written in terms of the pseudo-particle operators has 
the form: 



E 



dn k 
dt 



ijk 



ijk mjk 

+ Y T k c kcrlPma d J' J ~' J + E T k C kcr^ ~VW ~ h.C. 
mjk rnk 



(10) 



We set H = 1 and neglect changes in the electron spec- 
trum and local density of states in the tunneling contact 
leads, caused by the tunneling current. Therefore equa- 
tions of motion together with the constraint on the space 
of the possible system states (pseudo-particles number) 
(O give the following equations: 



ImJ2Tk-(c k ,f+b) 

ik 

ImY^Tk ■ (c ka df°-° fi- a ) 

ijk 

ImY,Tk-{ck„d+° a f la ) 

ijk 

Im E Tk ■ ( C kcr'>Pm- l 7 d j~ a ) 
mjk 



r fc^[(l - n ka (ei)) ■ nfia - n ka (£i) ■ n b ] 

i 

r fc / ,[(! _ n ka(Ejf j ' 7 ~ Ei-cr)) ■ n^ a - n ka (Ej Ij a - Ei-a) ■ 
ij 

Ik - n ka {EH } - E ia )) ■ n% - n ka {Ef^ - e la ) ■ n fi<T ] 

ij 

Tfe J^[(l — n ku (E" 1 [I cr — Ejjj' 7 )) ■ n^m-o- — n ka (E 7 / l II a — E° ' •■" 



J Uj i ' n dj 



mjk 



Im'Y T k ■ {Cka^t^ma) = T fe ^[(l - fl ka (E IV i - E^fj)) ■ n v - n ka (E Ivl - EjYj) ■ n^ma] 



(11) 



rnk 



Tunneling current I ka is determined by the sum of the 
right hand parts of the equations (fTTj) . 

Stationary system of equations can be obtained for the 



r 



pseudo particle filling numbers n/j, n^- a , , n^ m and 
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n = ^f. 
dt 



o = 



dn 



ipraa 



dt 



dn° 



dt 



= 



dt 



dun, 



fia 



Ot 



~T k Yy~n4, ma ■ n ka {Ei V i - Effi. - s k ) + n v • (1 - n ka (E Ivl - E^fi - e*))] + (k p) 

mcr 

-r fc • (1 - n^EJfj - E°- - e k )) - n^E^ - E°- - e k ) ■ n a -\ - 

3 

T k £[(1 - - Eg- - e k )) ■ n^ ma - ng- ■ n kt ,(Ey% - Eg- - e k )\ - 

3 

r*[-(l - n ka {Eivi - EY% - E k )) ■ n v + n^ ma ■ n ka (E IV i - Ejfj - e k )] + (k «-> p) 
-r fc ^[(1 - n fcCT (£^ - e< - E fc )) • </ - n ka (E™ - Ei - E k ) ■ n fla ] - 

i 

r k Y}n k - a (E™- - Ef% - E k ) ■ n% - (1 - n fe _ CT (^7}7 CT - Effi - s k )) ■ n^ m . a ] + (k o p) 
-r fc - n k - a {E°— - Ei ~ e k )) ■ n a — - n k . a (Ef— - Ei - e k ) ■ n fitr ] - 

iu 

T k J2i n kAEm - E a - - E k ) ■ rf— - (1 - n^iEfZ - E a - - e k )) ■ n^ a \ + (k o p) 

mtr 

rfe[n feCT (e l - e fe ) • rib - (1 - njtafo - ■ n /icr ] + 

T k £[(1 - n k ^(EJ— - Ei - E k )) ■ n a — - n k ^(EJ— - e 4 - e fc ) • n /m ] + 
T k - - et - e fe )) • - nhaiEfff -Ei- e k ) ■ n fia ] + (k o p) 

i 



(12) 



In these equations we neglect the non-diagonal aver- 
ages of pseudo-particle operators such as (f+bd + f a ) etc.. 
These terms are of the next order in small parameter 
T k ( p )/AE where AE is the energy difference between 
any energy states in the coupled QDs. We consider the 
paramagnetic situation, when conditions n/j CT = n/i_ CT , 
n^ma = n^m-v, n ka = n k - a and rig— = njg 7 are ful- 
filled. System of equations ([12")) in the stationary case 
is the linear system, which allows to determine pseudo 
particle filling numbers, electron filling numbers n e i(eV) 
and tunneling current I ka . 

It is necessary to mention that tunneling current 
through the proposed system can be also analyzed in 
usual terms of the electrons creation/annihilation opera- 
tors af a /ai a and c~l a /c ka in the localized and continuous 
spectrum states correspondingly: 



I = ha = ha = tk((c kcr a lrT } - (af a c krT )) 



(13) 



By means of Heisenberg equations of motion one can 
get system of equations exactly taking into account cor- 
relations of electron filling numbers in localized states in 
all orders [271, [28] (for weak tunneling coupling to the 
leads). But it is rather tedious to restore the informa- 
tion about the definite charge and spin configurations 
with different number of electrons from all order correla- 
tors for initial levels occupation numbers. So the method 



based on the pseudo particle operators is more convenient 
if we are interested in occupation of multi-electron states 
with particular charge and spin configuration. 



III. RESULTS AND DISCUSSION 

The behavior of the total electron occupation of the 
coupled QDs n e i (eV) and I — V characteristics strongly 
depends on the parameters of the tunneling contact: en- 
ergy levels position, the value of the Coulomb interaction 
and the relation between tunneling rates. The general 
features of the obtained results is the step- like I— V char- 
acteristics with non-equidistant steps related to the ener- 
gies of the different multi-electron states in the QDs and 
multiple charge redistribution between the two electron 
states with different spin configurations (singlet state and 
triplet state) which appears for the particular range of 
the system parameters and applied bias voltage. 

We first analyze the behavior of the the total electron 
occupation of the QDs n e i(eV) and I —V characteristics 
of the considered system for different single electron levels 
positions relative to the sample Fermi level and various 
tunneling rates to the contact leads (Fig. [T][3]). The bias 
voltage in our calculations is applied to the sample. Con- 
sequently, if both single electron levels are above(below) 
the Fermi level, all the specific features of the total elec- 
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FIG. 1: Total electron occupation of the QDs a).-c). and tunneling current d).-f). as a functions of the applied bias voltage 
for both single electron energy levels located above the sample Fermi level. Parameters e/Y = 0.8, T/T — 0.7, U/T = 1.85 are 
the same for all figures. a).,d). Y k = 0.01, Y p = 0.01; b).,e). Y k = 0.1, Y p = 0.01; c).,f). Y k = 0.01, Y p = 0.02. 



tron occupation and tunneling current characteristics can 
be observed at negative (positive) values of eV. In the 
case when both single electron energy levels are situated 
above the sample Fermi level (FigO} we observe the step- 
like behavior of the total electron occupation both for the 
symmetric (T^ — T p ) and asymmetric (r^ < T p ) tunnel- 
ing contact (Fig[Th,,c). The width and height of the steps 
are determined by the relation between the system pa- 
rameters T, e and U and r = Ffc + T p . 

In the particular range of the applied bias when con- 
dition T k > T. p is fulfilled the total electron occupation 
demonstrate significant jumps (FigUJa). 

The tunneling current is depicted in Figures [TJi-f as 
a function of the applied bias (tunneling current am- 
plitudes are normalized to 2 r fc + p ). In the asymmet- 
ric tunneling contact when condition < T p is valid 
(FigE) the tunneling current dependence on the ap- 
plied bias has a step-like structure. In the presence of 
Coulomb interaction negative tunneling conductivity ap- 
pears even for symmetric tunneling contact (FigJTJi). The 
negative differential conductivity is strongly pronounced 
in the asymmetric tunneling contact when the condition 
r fc > T p is fulfilled (FiglTk). 

In a QD without Coulomb interaction the total system 
occupation can only increase as applied bias increases and 
passes single electron levels. Quite different situation oc- 
curs in a system with strong Coulomb correlations. In 
the case when one of the single electron energy levels or 
both of them are situated below the sample Fermi level 
(Figl^JEJ) the total electron occupation of QDs demon- 



strates pronounced decreasing with increasing of the ap- 
plied bias (Fig. [2^-c; Fig[3^-c). Equilibrium occupation 
of two-electron states for zero bias leads to average occu- 
pation of one electron per spin (total occupation equal to 
2). But when the increasing bias reaches the energies of 
multi-electron excited states, the total occupation begins 
to decrease. Using single-electron language we can say 
that additional tunneling electrons "push out" electrons 
from the states below the Fermi level due to Coulomb 
repulsion. Or one can look at this effect as increasing 
of probability for electrons to leave the QD due to ap- 
pearance of several non-elastic channels of tunneling (ac- 
companied with changing of multi-electron states of the 
QD). 

The tunneling current as a function of the applied bias 
for this case is depicted in Figures f OJi- f and reveals 
the monotonic step-like behavior. In both cases the upper 
single electron energy level (eq + T) does not appear as 
a step in the I — V characteristics. 

Another interesting feature which appears due to cor- 
relations is multiple charge redistribution between the 
two electron states with different spin configurations (sin- 
glet state and triplet state) as a function of the applied 
bias voltage (Figf4][6]). 

In the case of both single electron energy levels situ- 
ated above the sample Fermi level (Fig@| several ranges 
of applied bias exist where the triplet state occupation 
exceeds the singlet state occupation both for the strong 
{U/Y > T/T) and weak (U/Y ~ T/T) Coulomb interac- 
tion. In the case of strong Coulomb interaction localized 
charge is mostly accumulated in the triplet state. 




FIG. 2: The same dependencies as in the Fig[T]but one single electron energy level is located above and the other - below the 
sample Fermi level. Parameters s/T — —0.5, T/T — 0.8, U/T = 1.5 are the same for all figures. a).,d). Tt = 0.01, T p — 0.01; 
b).,e). F k = 0.02, F p = 0.01; c).,f). T k = 0.01, T p = 0.02. 



The occupation of the triplet state with the fixed value 
of the spin projection Sz is lower than the the occupa- 
tion of the singlet state if Coulomb repulsion is weak 
(U/T ~ T/T) (Fig|la-c). In the case of strong Coulomb 
interaction (U/T ^ T/T) four ranges exist where charge 
is quite equally distributed between the singlet state and 
triplet state with the fixed value of the spin projection 
(Fig|4ji,e,f) for different ratios between the tunneling 
transfer rates. 

Figure [5] demonstrates calculation results in the case 
when one single electron energy level is located above 
and the other - below the sample Fermi level. For the 
weak Coulomb interaction (U/T ~ T/r)the charge is 
quite completely localized in the singlet state both in the 
symmetric and asymmetric tunneling contact (Fig(5^-c). 

In the case of strong Coulomb interaction (U/T 3> 
T/T) charge in the system is also mostly located on the 
singlet state for the wide range of applied bias (FigJSJl-f), 
but for some ranges of the applied bias occupation of the 
triplet state exceeds occupation of the singlet state. But 
for symmetric tunneling contact the triplet state with 
the fixed spin projection is always less occupied than the 
singlet state . 

The situation when both single electron energy levels 
are positioned below the sample Fermi level is demon- 
strated in Fig|51 In the case of weak Coulomb interaction 
only one range of the applied bias exists where occupa- 
tion of the triplet state is equal to or even exceeds the 
occupation of the singlet state (FigJSJi-c). Occupation 
of the triplet state with the fixed spin projection is al- 
ways lower than the occupation of the singlet state. The 



increasing of the Coulomb interaction leads to the for- 
mation of several ranges of the applied bias where triplet 
state filling numbers exceed singlet state filling numbers 

(FigEH-f). 



IV. CONCLUSION 

We investigated tunneling through the system of two 
interacting QDs weakly coupled to the electrodes with 
Coulomb interaction between localized electrons. In the 
considered system if electrons number N is changed due 
to the tunneling processes the modification of the energy 
spectrum is not reduced to the simple adding of Coulomb 
interaction U per electron. One, two, three or four elec- 
trons can be localized in the coupled QDs , each state 
with fixed total charge and spin projection has it's own 
energy. Transitions between these states were analyzed in 
terms of pseudo-particle operators with constraint on the 
possible physical states of the system. Filling numbers 
of different multi-electron states, total electron occupa- 
tion of QDs and I — V characteristics were investigated 
for different single electron levels positions relative to the 
sample Fermi level and various tunneling transfer rates. 

It was shown that total electron occupation demon- 
strates in some cases significant decreasing with increas- 
ing of applied bias - contrary to the situation with no 
correlations. 

We revealed that for some parameter range, the system 
demonstrates negative tunneling conductivity in certain 
ranges of the applied bias voltage due to the Coulomb 
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FIG. 3: The dependencies of total electron occupation of the QDs a).-c). and tunneling current d).-f). on the applied bias 
voltage for both single electron energy levels located below the sample Fermi level. Parameters e/T = —0.5, T/T — 0.3, 
U/T = 1.0 are the same for all cases. a).,d). F k = 0.01, T p = 0.01; b).,e). F fc = 0.02, F p = 0.01; c).,f). F k = 0.01, T p = 0.02. 



correlations. A negative tunneling conductivity is well 
pronounced if both energy levels are located above the 
Fermi level. When energy levels are located on the op- 
posite sites of the Fermi level or both of them are posi- 
tioned below the Fermi level negative tunneling conduc- 
tivity was not observed. 

Coulomb correlations result in multiple charge redistri- 
bution between the two-electron states with different spin 
configurations (singlet and triplet states) with changing 
of applied bias voltage. It was found that for particular 



range of the system parameters the triplet-state occupa- 
tion can exceed the singlet-state occupation (the inverse 
occupation takes place). 

So tunneling properties of correlated electron systems 
can be correctly described only in terms of multi-electron 
states, which allows to find some unexpected effects. 

This work was partly supported by the RFBR and 
Leading Scientific School grants. The support from the 
Ministry of Science and Education is also acknowledged. 
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FIG. 4: Occupation of two electron state for different spin configurations as a functions of the applied bias voltage for both 
single electron levels located above the sample Fermi level. Filling numbers in the singlet state are shown by the black line, 
filling numbers in the triplet state with the fixed projection of the spin are shown by the grey line, full filling numbers in 
the triplet state are shown by the black-dashed line. Parameters e/Y = 0.5, T/Y = 0.3 are the same for all figures. a).,d). 
T fe = O.Of, Y p = O.Of; b).,e). Y k = 0.02, Y p = O.Of; c).,f). Y k = O.Of, Y p = 0.02. a).-c). U/Y = 0.5; d).-f). U/Y = 2.0. 




FIG. 5: The dependence of two electron state filling numbers for different spin configurations on applied bias voltage when one 
of the single electron energy levels is located above and the other- below the sample Fermi level. Filling numbers in the singlet 
state are shown by the black line, filling numbers in the triplet state with the fixed projection of the spin are shown by the grey 
line, full filling numbers in the triplet state are shown by the black-dashed line. Parameters e/Y = —0.3, T/Y = 0.5 are the 
same for all figures. a).,d). Y k = 0.01, Y p = 0.01; b).,e). Y k = 0.02, Y p = 0.01; c).,f). Y k = 0.01, Y p = 0.02. a).-c). U/Y = 0.5; 
d).-f). U/Y = 2.0. 
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FIG. 6: The same dependencies as on the Fig |4l5l but for both single electron energy levels located below the sample Fermi level. 
Filling numbers in the singlet state are shown by the black line, filling numbers in the triplet state with the fixed projection 
of the spin are shown by the grey line, full filling numbers in the triplet state are shown by the black-dashed line. Parameters 
e/r = -0.5, T/r = 0.3 are the same for all figures. a).,d). T k = 0.01, F p = 0.01; b).,e). F k = 0.02, T p = 0.01; c).,f). T k = 0.01, 
r p = 0.02. a).-c). U/r = 0.5; d).-f). U/Y = 2.0. 



